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ABSTRACT 

We present upper-limits on the masses of the putative central intermediate-mass black 
holes in two nearby Galactic globular clusters: 47Tuc (NGC104), the second brightest 
Galactic globular cluster, and NGC6397, a core-collapse globular cluster and, with 
a distance of 2.7 kpc, quite possibly the nearest globular cluster, using a technique 
suggested by T. Maccarone. These mass estimates have been derived from 3cr upper 
limits on the radio continuum flux at f.4 GHz, assuming that the putative central 
black hole accretes the surrounding matter at a rate between 0.1 % and 1 % of the 
Bondi accretion rate. For 47Tuc, we find a 3tr upper limit of 2060 — 670M Q , depending 
on the actual accretion rate of the black hole and the distance to 47Tuc. For NGC6397, 
which is closer to us, we derive a 3cr upper limit of 1290 — 390M Q . While estimating 
mass upper-limits based on radio continuum observations requires making assumptions 
about the gas density and the accretion rate of the black hole, their derivation does not 
require complex and time consuming dynamical modeling. Thus, this method offers 
an independent way of estimating black hole masses in nearby globular clusters. If, 
generally, central black holes in stellar systems accrete matter faster than 0.1 % of the 
Bondi accretion rate, then these results indicate the absence of black holes in these 
globular clusters with masses as predicted by the extrapolated M. — a c relation. 
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1 INTRODUCTION 

Theory provides us with different ways of producing 
intermediate-mass black holes (IMBHs) in dense stellar clus- 
ters, such as globular clusters. IMBHs are defined here as 
having a mass larger than that of stellar black holes but 
smaller than that of the super-massive black holes (SMBH) 
found in galaxies, so 10 2 Mq < M, < 10 6 Mq. One scenario 
envisages an IMBH to grow from a stellar-mass seed black 
hole, the remnant of a massive star that exploded as a super- 
nova (Miller & Hamilton 2002). This seed black hole rapidly 
sinks to the cluster center due to dynamical friction, or, in 
the case of a post-core-collapse cluster, was formed near the 
cluster center since that is where the most massive stars can 
be found, and grows by accreting stars, less massive stellar 
black holes or interstellar gas. In a core-collapse cluster, stars 
near the core are packed so densely that some may collide, 
causing a runaway growth of the remnant's mass (Portc- 
gies Zwart & McMillan 2002). Both mechanisms seem to 
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lead naturally to IMBHs with a mass that is about 0.1% 
of the cluster mass. The IMBH can be even more massive 
(up to 1O 4 M0) if it originated from multiple seeds. Another 
suggested way of forming IMBHs considers the radiation 
drag that is exerted by stars on the interstellar medium 
(Kawakatu et al. 2003; Kawakatu & Umemura 2005). Ra- 
diation drag is expected to funnel gas, expelled by SNu and 
SNla explosions, towards the cluster center where it is likely 
to form a massive central object surrounded by an accre- 
tion disk. It is expected that this central object will collapse 
into an IMBH if its mass exceeds a few 100 Mq. In massive 
galaxies, this process may lead to central black holes that 
have a mass of about 0.1% of the bulge mass (Kawakatu et 
al. 2003). However, in the case of globular clusters, IMBHs 
formed this way are expected to be much less massive than 
0.1% of the cluster mass (the ratio of the IMBH mass to the 
cluster mass is rather of the order of 10 _J ). In fact, clusters 
less massive than ~ 5 x 1O 6 M0 are not expected to contain a 
central black hole, although they can harbor a central object 
not massive enough to collapse into a black hole (Kawakatu 
& Umemura 2005). Yet another, somewhat more specula- 
tive, possibility is that Population ill stars, with masses of a 
few hundred solar masses, were the seeds of IMBHs (Madau 
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& Rocs 2001). Nuclear energy production in Population ill 
stars more massive than ~ 260Mq cannot counteract grav- 
ity and such objects immediately collapse into black holes. 
Such very massive seed black holes, which could be called 
IMBHs in their own right, then grow further by accretion. 
This way of producing IMBHs, unlike the previous ones, 
does not require the host stellar system to be dense. 

Therefore, even upper limits on the masses of central 
black holes in globular clusters are of great interest, since 
they constrain the possible formation scenarios for IMBHs 
and the nature of IMBH progenitors. In this paper, we 
present upper limits on the IMBH masses of 47Tuc and 
NGC6397, two of the nearest globular clusters, derived from 
20 cm continuum observations, obtained with the Australia 
Telescope Compact Array (ATCA). In Section 2, we briefly 
discuss how IMBHs can be detected via their radio con- 
tinuum emission. In Sect. 3, the observations and the data 
reduction procedures are presented, followed by a discussion 
of the results in Sect. 4. We summarize our conclusions in 
section 5. 



2 DETECTING IMBHS USING RADIO 
CONTINUUM EMISSION 

While various approaches to the estimation of IMBH masses 
in globular clusters have been suggested and/or applied 
(such as dynamical modeling of stellar kinematics, see e.g. 
Gerssen et al. (2002), or the detection of high-proper- motion 
stars near the globular cluster center (Drukier & Bailyn 
2003)) one of the observationally least demanding methods 
is the possibility to detect black holes through the radio 
emission they produce. Gas that is accreted onto a black 
hole converts part of its rest mass into radiation which, in 
the mass regime of IMBHs, is most easily recovered at radio 
wavelengths. The equation that relates the mass of a black 
hole, M,, in a stellar system at a distance d to its X-ray 
luminosity Lx and radio continuum emission F con t reads as 



= 10 



3 x 10 31 ergs" 1 
d 



M. 



100 M, 



lOkpc 



(1) 



assuming a flat radio spectrum (Maccarone 2004) and the 
Fundamental-Plane relationship between radio luminosity, 
X-ray luminosity, and black hole mass (Merloni, Heinz, di 
Matteo 2003; Falcke, Kording, Markoff 2004). For a given 
accretion rate, which is usually expressed as a fraction of 
the Bondi accretion rate 
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(with n and T the gas density and temperature, respectively) 
the X-ray luminosity can be calculated, assuming that 10 % 
of the rest mass energy of the infalling matter is converted 
into radiation, and that the total luminosity L to t is related 
to the X-ray luminosity as given by Ltot = 6 x 10 _3 Lx 8 + 
Lx (using the Eddington luminosity as unit). This relation 
takes into account the energy of a relativistic jet (Fender 



et al. 2003). This way, the radio continuum flux F con t can 
be expressed as a function of the black hole mass M,, the 
distance d, and the gas density n and temperature T. In the 
following, we will assume n = 0.067 cm" " 3 and T — 10 4 K 
(see Freire et al. (2001)) as typical values. 

As shown by e.g. Di Matteo et al. (2001) and Quataert 
& Gruzinov (2000) for nearby early-type galaxies and for the 
Milky Way, respectively, massive central black holes seem to 
accrete much slower than predicted by the Bondi estimate 
(with M. < 10~ 2 M., Bondi)- The low number of Galactic 
neutron stars that were detected in the ROSAT All-Sky 
Survey can be explained if they are accreting at a rate of 
roughly 0.1 % of the Bondi rate (Perna et al. 2003), al- 
though, of course, a neutron star accreting from the inter- 
stellar medium is a far cry from a central black hole. Theo- 
retical models of the accretion of matter that has non-zero 
angular momentum onto a black hole (Krumholz, McKee, 
Klein 2005; Proga & Begelman 2003) show a reduction of 
the accretion rate down to the level of a few percent of the 
Bondi rate. In the following, we will assume an accretion 
rate in the range 10~ 3 M., Bond i < M. < 10~ 2 M., Bond i, in 
agreement with Maccarone et al. (2004) and Maccarone et 
al. (2005). Thus, it is possible to estimate the masses (or 
place upper limits thereon) of central black holes in nearby 
globular clusters using straightforward and relatively short 
radio observations. 



3 OBSERVATIONS AND DATA REDUCTION 

We observed both NGC6397 and 47Tuc with the ATCA on 
respectively 2004 December 23 and 26. We used the 1.5D 
configuration, with baselines ranging from 107 to 4439 m. 
The continuum observations were made at a central fre- 
quency of 1.384 GHz in the FULL_128_1 correlator con- 
figuration with a bandwidth of 128 MHz divided over 32 
channels into 4 products. At the start of each observation, 
we observed the standard ATCA primary calibrator 1934- 
638 for 10 min. A secondary calibrator was observed ev- 
ery 45 min. The total integration time (including the cal- 
ibration) for each globular cluster was 3.5 h. The data re- 
duction was performed with the MIRIAD package (Sault, 
Teuben, Wright 1995), the standard ATCA data analysis 
program. The observed data were loaded into MIRIAD with 
the birdie, xycorr, and reweight options, which respectively 
flags out the channels that are affected by the ATCA self- 
interference, corrects for the phase difference between the 
X and Y channels and re-weights the visibility spectrum in 
the lag domain to eliminate the effects of the Gibbs phe- 
nomena. This results in 13 channels, each 4 MHz wide. The 
usual data reduction steps, including phase, amplitude and 
bandpass calibration, were afterwards performed. Instead of 
applying the traditional CLEAN algorithm (Conway, Corn- 
well, Wilkinson 1990) we used the Multi-Frequency Synthe- 
sis (MFS) method (Sault & Conway 1999) with a natural 
weighting to surpress the noise. This improved drastically 
the (u, w)-coverage by eliminating spectral artifacts. 

For NGC6397, after applying all above reduction steps, 
we constructed a radio continuum map with a beam size of 
36.08" x 15.64" and a la noise value of 72 /ijy. No radio 
continuum sources could be detected in this object at this 
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noise level. We therefore assume 216 /xjy as the 3<r upper 
limit on the radio emission of the black hole. 

Our 20 cm map of 47Tuc has a depth similar to that 
presented by McConnell et al. (2004) although it has higher 
spatial resolution : our final image has a beam size of 
33.39" x 12.55" and a la noise value of 75 pJy. None of the 
radio continuum sources is close enough to the cluster cen- 
ter to be associated with a putative black hole. We therefore 
assume 225 /iJy as the 3a upper limit on the radio emission 
of the black hole. 



4 RESULTS 

47Tuc is a massive globular cluster, with a total stellar mass 
M » 1.5 x 1O 6 M (Gebhardt & Fischer 1995). Its large core 
radius suggests that it is a relaxed (and not a post-core- 
collapse) cluster (Howell, Guhathakurta, Gilliland 2000) . We 
use the distance modulus m — M = 13.50 ± 0.08 mag de- 
rived by Gratton et al. (2003). NGC6397 on the other hand 
is a much less massive cluster, with a total stellar mass 
M » 1.0 x 10 5 M Q (Meylan & Mayor 1991). It is a very cen- 
trally concentrated, collapsed cluster. Using accurate color- 
magnitude diagrams and luminosity functions of two fields 
inside NGC6397, Andreuzzi et al. (2004) found strong ev- 
idence for different mass distributions and hence for mass 
segregation to have occurred. We adopt a distance modulus 
m-M = 12.13±0.15 mag (Reid & Gizis 1998), determined 
by main-sequence fitting using lower main sequence stars 
with hipparcos parallaxes measured to a precision of better 
than 10%. 

We estimate upper-limits on the masses of the putative 
central IMBHs in these globular clusters from 3a upper lim- 
its on the radio continuum flux. We assume that the putative 
central black hole accretes the surrounding matter at a rate 
between 0.1 % and 1 % of the Bondi accretion rate and take 
into account the uncertainty on the distance modulus. For 
47Tuc, we find a 3er upper limit for the mass of the central 
black hole of 2060 - 67OM . For NGC6397, we derive a 3a 
upper limit of 1290 — 390M© for the black hole mass. 

Massive elliptical and spiral galaxies seem to adhere 
closely to a relation between the mass of the central black 
hole, M. , and the central velocity dispersion, a c (see Fig. 1). 
Depending on the sample selection and the way the linear 
regression is done, the slope and zero-point of this relation 
vary (Ferrarese & Ford 2005; Tremaine, Gebhardt, Bender 
2002). The existence of this relation is however generally 
accepted. Although there is theoretically no obvious rea- 
son why this relation should hold with the same slope also 
in the mass regime of globular clusters, with the growth 
of SMBHs most likely being driven by the merger histo- 
ries of their host galaxies, it is nevertheless tempting to do 
so. The central velocity dispersion of 47Tuc, ao ~ 12 km/s 
(Gebhardt & Fischer 1995), corresponds to an IMBH mass 
M. = 1600 M Q . The central velocity dispersion of NGC6397 
is very low, ao ~ 5.0 km/s (Meylan & Mayor 1991), and cor- 
responds to an IMBH mass M m — 50 Mq . 

The upper limit for the mass of the black hole in 47Tuc 
lies just below the extrapolated M. — a c relation of mas- 
sive elliptical and spiral galaxies as derived by Tremaine et 
al. (2002). As this concerns a 3a upper limit, this means 
that the mass of the central black hole in 47Tuc must be 
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Figure 1. The M. — <r c relation of massive elliptical and spi- 
ral galaxies, dwarf galaxies, and globular clusters. The M. — a c 
derived by Ferrarese & Ford (2005) lies within the red region, 
which quantifies the la uncertainty on that relation. The rela- 
tion and its la uncertainty as derived by Tremaine et al. (2002) 
is indicated in green. The blue datapoints correspond to the 
3a upper limits on the masses of the black holes in 47Tuc and 
NGC6397, the upper dot corresponding to an accretion rate equal 
to M. = 10~ 3 M. i B ondi , the lower to M. = 10 _2 M., B ondi (the 
same holds for ajCen). 



much lower than predicted by this relation, on the condi- 
tion, of course, that the black hole accretes at a rate higher 
than 0.1 % of the Bondi rate. It also excluded the pres- 
ence of a IMBH with a mass of about 0.1 % of the clus- 
ter mass. The extrapolated M. — a c relation as derived by 
Ferrarese & Ford (2005) predicts lower M.-values in the 
globular-cluster regime than the Tremaine et al. (2002) ver- 
sion and is potentially consistent with our results. The up- 
per limit for the mass of the black hole in NGC6397 lies 
above both renditions of the M, — a c relation and thus puts 
no extra constraints on this relation in the globular-cluster 
regime. Given the uncertainty concerning the detection of 
a black hole in M15 (Gerrsen et al. 2002; McNamara, Har- 
rison, Anderson 2003; Maccarone, Fender, Tzioumis 2005), 
this means that no secure detection of an IMBH exists in 
the low-mass regime (a c < 15 km/s) of the M. — a c relation. 
Radio-continuum observations of uiCen yielded a 3a upper 
limit of M, = 1470 — 460, depending on the accretion rate 
(Maccarone, Fender, Tzioumis 2005). We adopt a distance 
of 4.8±0.3 kpc and a central velocity dispersion 18±2 km/s 
(van de Ven et al. 2006). Dynamical modeling of the central 
regions of the Local Group galaxies NGC205 (Valluri et al. 
2005) and M33 (Merritt, Ferrarese, Joseph 2001; Gebhardt 
et al. 2001) has yielded no conclusive evidence for IMBHs 
in the a c ~ 30 km/s regime (with M33, in fact, well below 
the la uncertainties of the extrapolated M. — a c relations) . 
This leaves Gl as the only object with a relatively secure de- 
tection of a sub-10 6 M IMBH (Gebhardt, Rich, Ho. 2002; 
Baumgardt et al 2003; Gebhardt, Rich, Ho 2005) and that 
agrees with the extrapolated M. — a c relation. 
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5 CONCLUSIONS 

Using radio-continuum observations with the ATCA radio 
telescope array, we estimate the 3a upper-limits on the 
masses of the putative central IMBHs in two nearby Galac- 
tic globular clusters, 47Tuc and NGC6397, at respectively 
2060 - 670M© and 1290 - 390M Q . These mass estimates 
have been derived from 3a upper limits on the radio contin- 
uum flux. We assume that the putative central black hole 
accretes the surrounding matter at a rate between 0.1 % and 
1 % of the Bondi accretion rate and take into account the 
uncertainty on the distance modulus. 

Black hole masses estimated using radio-continuum ob- 
servations by necessity depend on assumptions regarding 
the density and temperature of the interstellar medium, and 
on the accretion rate. Observations of a few globular clus- 
ters provide "typical" values for the characteristics of the 
interstellar medium. The black hole accretion rate can be 
constrained by observations of the central black hole of the 
Milky Way and by theoretical models of accretion disks. 
Still, radio-continuum estimates of black hole masses prob- 
ably cannot boast the same degree of accuracy as mass esti- 
mates from e.g. dynamical models. However, while dynami- 
cal models fitted to radial velocities and proper motions of 
stars near the center of a globular cluster yield the most 
accurate mass estimates, their spatial resolution is not high 
enough to identify the central dark mass as a black hole. If, 
on the other hand, one detects radio emission coming from 
accreted matter falling into a steep gravitational well then 
one can be sure that the central object must be very com- 
pact, i.e. a black hole. A cluster of white dwarfs or neutron 
stars would not be able to generate such a feature. In a way, 
both methods are complementary. Still, it seems telling that 
we obtain a 3a upper-limit for the mass of the putative cen- 
tral black hole in 47Tuc, a nearby massive globular cluster, 
that places this object marginally below the extrapolation 
of the M. — a c relation of the bright ellipticals and spirals. 
The data for NGC6397 unfortunately do not constrain the 
low-mass end of the M. — a c relation. 

Our results hopefully form the basis for further at- 
tempts to detect IMBHs in these globular clusters, which, in 
the case of a detection, would lead to a better understanding 
of the low-mass tail of the black-hole population and of the 
accretion rates of central IMBHs in globular clusters. 
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